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Abstract 

This paper is devoted to the study of the performance of the Linear Minimum Mean-Square Error receiver for 
(receive) correlated Multiple-Input Multiple-Output systems. By the random matrix theory, it is well-known that the 
Signal-to-Noise Ratio (SNR) at the output of this receiver behaves asymptotically like a Gaussian random variable as 
the number of receive and transmit antennas converge to +00 at the same rate. However, this approximation being 
inaccurate for the estimation of some performance metrics such as the Bit Error Rate and the outage probability, 
especially for small system dimensions, Li et al. proposed convincingly to assume that the SNR follows a generalized 
Gamma distribution which parameters are tuned by computing the first three asymptotic moments of the SNR. In this 
article, this technique is generalized to (receive) coiTelated channels, and closed-form expressions for the first three 
asymptotic moments of the SNR are provided. To obtain these results, a random matrix theory technique adapted to 
matrices with Gaussian elements is used. This technique is believed to be simple, efficient, and of broad interest in 
wireless communications. Simulations are provided, and show that the proposed technique yields in general a good 
accuracy, even for small system dimensions. 

Index Terms: Large random matrices, correlated channels, outage probability. Bit Error Rate (BER), Gamma 
approximation, minimum mean square error, Multiple-Input Multiple-Output (MIMO) systems, Signal-to-Noise 
Ratio (SNR). 

I. Introduction 

Since the mid-nineties, digital communications over Multiple Input Multiple Output (MIMO) wireless channels 
have aroused an intense research effort. It is indeed well-known since Telatar's work [1] that antenna diversity 
increases significantly the Shannon mutual information of a wireless link; In rich scattering environments, this 
mutual information increases linearly with the minimum number of transmit and receive antennas. Since the findings 
of [1], a major effort has been devoted to analyse the statistics of the mutual information. Such an analysis has 
strong practical impacts; For instance, it can provide information about the gain obtained from scheduling strategies 
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[2]; it can be used as a performance metric to optimally select the active transmit antennas [3], etc. 
The early results on MIMO channels mutual information concerned channels with centered independent and 
identically distributed entries. It is of interest to study the statistics of this mutual information for more practical 
(correlated) MIMO channels. In this course, many works established the asymptotic normality of the mutual 
information in the large dimension regime for the so called Kronecker correlated channels [4], [5], for general 
spatially correlated channels [6] and for general variance profile channels [7]. 

Another performance index of clear interest is the Signal to Noise Ratio (SNR) at the output of a given receiver. 
In this paper we focus on one of the most popular receivers, namely the linear Wiener receiver, also called LMMSE 
for Linear Minimum Mean Squared Error receiver. In this context, an outage event occurs when the SNR at the 
LMMSE output lies beneath a given threshold. One purpose of this paper is to approximate the associated outage 
probability for an important class of MIMO channel models. Another performance index associated with the SNR 
is the Bit Error Rate (BER) which will be also studied herein. 

Outage probability approximations has been provided in recent works for various channels, under very specific 
technical conditions (in the case where the moment generating function [8] or the probability density function [9| 
have closed form expressions; when a first order expansion of the probabihty density function can be derived [lOJ; 
in the more general case where the moment generating function can be approximated by using Pade approximations 
[1 1]; etc.). All these results deal with specific situations where the statistics of the SNR could be derived for finite 
system dimensions. 

Alternatively, by making use of large random matrix theory, one can study the behavior of the SNR in the asymptotic 
regime where the channel matrix dimensions grow to infinity. For fairly general channel statistical models, it is then 
possible to prove the convergence of the SNR to deterministic values and even establish its asymptotic normality 
(see for instance [12], [13]). However, this Gaussian approximation is not accurate when the channel dimensions 
are small. This is confirmed in e.g. [14] where it is shown that the asymptotic BER based on the sole Gaussian 
approximation is significantly smaller than the empirical estimate. A more precise approximation of the BER or the 
outage probability is expected if one chooses to approximate the SNR probability distribution with a distribution 
1) which is supported by M-|_ (indeed, a Gaussian random variable takes negative values which is not realistic), 2) 
which is adjusted to the first three moments of the SNR instead of the first two moments needed by the Gaussian 
approximation. 

In this line of thought, Li, Paul, Narasimhan and Cioffi [15] proposed to use alternative parameterized distributions 
(Gamma and generalized Gamma distributions) whose parameters are set to coincide with the asymptotic moments 
of the output SNR. This approach was derived for (transmit) correlated channels and asymptotic moments were 
provided for the special case of uncorrelated or equicorrelated channels. For the general correlated channel case, 
only limiting upper bounds for the first three asymptotic moments were provided. Based on Random Matrix Theory 
and especially on the Gaussian mathematical tools elaborated in [4] and further used in [16], we derive closed-form 
expressions for the first three moments, generalizing the work of [15] to a general (receive) correlated channel. 
Using the generalized Gamma approximation, we provide closed-form expressions for the BER and numerical 



October 16, 2008 



DRAFT 



3 



approximations for the outage probability. 
Paper organization 

In section II, we present the system model and derive the SNR expression. Then we review in section III the 
Generalized Gamma approximation before providing the asymptotic central moments in the next section. Finally, 
we discuss in the last section the simulation results. 

II. System Model and SNR expression 

We consider an uplink transmission system, in which a base station equipped by N correlated antermas detects 
the symbols of a given user of interest in the presence of K interfering users. The N dimensional received signal 
writes: 

r = Ss + n, 

T 

where s = [sq, • • • , sr] is the transmitted complex vector signal with size K + 1 satisfying Ess* = Ik+i, and 
S is the A?^ X {K + 1) channel matrix. We assume that this matrix writes as 

where ^ aNxN Hermitian nonnegative matrix that captures the correlations at the receiver, P = diag (po > • ■ • jPk) 
is the deterministic matrix of the powers allocated to the different users and W = [wo, • • • , w^] (wfc being the 
fcth column) is a A'^ x {K + 1) complex Gaussian matrix with centered unit variance (standard) independent and 
identically distributed (i.i.d) entries. To detect symbol sq and to mitigate the interference caused by users 1,. . . ,K, 
the base station applies the LMMSE estimator, which minimizes the following metric: 

g = minE |h*r — SqI^ • 

h 

Let y = then it is well known that the LMMSE estimator is given by: 

g= (EE* +pljv)"'y- 

Writing the received vector r = soy + where soY is the relevant term and represents the interference plus 
noise term, the SNR at the output of the LMMSE estimator is given by : Pk = |g*y|^ /IE |g*rin|^. Plugging the 
expression of g given above into this expression, one can show that the SNR is given by: 

1 



= y* WPW**5 + pl^j y, 

with P = diag (pi, • • • ,Pk) and W = [wi, • • • , wk]- Let ^ = UDU* be a spectral decomposition of Then, 
(3k writes: 

(3k = ^wjUD^ ^-^D5U*WPW*UD5 +pl;v^ D^U^wo , 
pK \Kp 
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where: z = U*wo (resp. Z = U*W) is a iV x 1 vector with complex independent standard Gaussian entries (resp. 
N X K matrix with independent Gaussian entries). 

Under appropriate assumptions, it can be proved that admits a deterministic approximation as K,N ^ oo, 
the ratio being bounded below by a positive constant and above by a finite constant. Furthermore, its fluctuations 
can be precisely described under the same asymptotic regime (for a full and rigorous computation based on random 
matrix theory, see [13]). As it will appear shortly, a deterministic approximation of the third centered moment of 
I3k is needed and will be computed in the sequel. 

III. Bit Error Rate and outage probability approximations 

A. A quick reminder of the generalised Gamma distribution 

Recall that if a random variable X follows a generalized gamma distribution G{a,b,£), where a and h are 
respectively referred to as the shape and scale parameters, then: 



The probability density function (pdf) of the generalized Gamma distribution with parameters {a, b, Q does not 
have a closed form expression but its moment generating function (MGF) writes: 



B. BER approximation 

Under QPSK constellations with Gray encoding and assuming that the noise at the LMMSE output is Gaussian, 
the BER is given by: 



where (3^ denotes an asymptotic deterministic approximation of the first moment of Pk- It was shown however in 
[15] that this expression is inaccurate since a Gaussian random variable aUows negative values and has a zero third 
moment while the output SNR is always positive and has a non-zero third moment for finite system dimensions. 
To overcome these difficulties, Li et al. [15] approximate the BER by considering first that the SNR foUows a 
Gamma distribution with scale a and shape 6, these parameters being tuned by equating the first two moments of 
the Gamma distribution with the first two asymptotic moments of the SNR. However, the third asymptotic moment 
was shown to be different from the third moment of the Gamma distribution which only depends on the scale a and 
shape b. In light of this consideration, Li et al. [15] refine this approximation and consider that the SNR foUows 
a generalized Gamma distribution which is adjusted by assuming that its first three moments equate the first three 



EX = ab, var(X) = ab'^ and E{X - EXf = + l)ab^ . 





) 
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asymptotic moments of the SNR. As expected, this approximation has proved to be more accurate than the Gamma 
approximation, and so will be the one considered in this paper Next, we briefly review this technique, which we 
will rely on to provide accurate approximations for the BER and outage probability. 

Let Eoc{Pk), varoo(/3/f) and Soo{Pk) denote respectively the deterministic approximations of the asymptotic 
central moments of Pk- Then, the parameters ^, a and b are determined by solving: 

^oo{Pk) = ab, vaiooiPK) = ab"^ and Soo(/?k) = (C + 

thus giving the following values: 

^ _ {E^if3K)? ^_ vareo(/3K) ^ _ Soo(fe)E^(fe) ^ 

varoo(/3K) ' Eoo(/3k) (varoo(/3A:))^ 

Using the MGF, one can evaluate the BER by using the following relation [19], that holds for QPSK constellation: 

TT Jo V 2 sin^ 

Note that similar expressions for the BER exist for other constellations and can be derived by plugging the following 
identity involving the function Q{x) [19]: 

TT Jo \ 2 am 0/ 

into the BER expression. 

C. Outage probability approximation 

Only the moment generation function (MGF) has a closed form expression. Knowing the MGF, one can compute 
numerically the cumulative distribution function by applying the saddle point approximation technique [20]. Denote 
by K{y) = log(MGF(y)) the cumulative generating function, by y the threshold SNR and by ty the solution of 



BER = - / MGF ( . ^ , ) d(j>. (1) 



K'{ty) = y. Let wo and uq be given by: wo = sign{ty) ^y 2 [tyy — K{ty)) and uq — tyy/JCUlyj. The saddle point 
approximate of the outage probability is given by: 

Pout = $K) + Hwo) (—-—], (2) 



,m Uo, 

where $(a;) = vl7^~* '^^ ^"'^ "^(^^ ~ V^^^^ denote respectively the standard normal cumulative 
distribution function and probability distribution function. 

So far, we have presented the technique that will be used in simulations for the evaluation of the BER and outage 
probability. This technique is heavily based on the computation of the three first asymptotic moments of the SNR 
Pk, an issue that is handled in the next section. 

IV. Asymptotic moments 

A. Assumptions 

Recall from Section the various definitions ivT, TV, D,D. In the following, we assume that both K and N go 
to +00, their ratio being bounded below and above as follows: 

K K 
< ^~ = liminf — < £+=limsup — < +oo . 
N ~ N 
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In the sequel, the notation K ^ co will refer to this asymptotic regime. We will frequently write "Dk and Dx to 
emphasize the dependence in K, but may drop the subscript K as well. Assume the following mild conditions: 
Assumption Al.- There exist real numbers dmax < oo and dmax < oo such that: 

sup IIDkII < dmax and sup IjDxIl < dmax, 

K K 

where ||Dx|| and ||Di^|| are the spectral norms of Dx and Dx- 
Assumption A2: The normalized traces of D^f and D^e satisfy: 

inf 4Tr(DK) > and inf ^TrfD/^) > 0. 

K K K K 

B. Asymptotic moments computation 

In this section, we provide closed form expressions for the first three asymptotic moments. We shall first introduce 
some deterministic quantities that are used for the computation of the first, second and third asymptotic moments. 
Proposition 1: (cf. [4]) For every integer K and any < > 0, the system of equations in ((5, 5) 

Ik = ^TrD,^ (l + tJ^D^' 

admits a unique solution [^SK{t),5K{t)^ satisfying (5/f(t) > 0, ^x(t) > 0. 
Let T and T be the N x N and K x K diagonal matrices defined by: 



T = ( I + tSKD ] ^ and T = ( I + tSK^ 



Note that in particular: S = ;^TrDT and S = -^TrDT. Define also 7 and 7 as 7 = -^TrD^T^ and 7 = -^TrD^T^. 
Finally, replace t by ^ and introduce the following deterministic quantities: 

o2 _ 7 f 11 



Vk 



P -11 

2p3 



~\3 



p3 



As usual, the notation uk = 0{(3k) means that aKif^K)^^ is uniformly bounded as X ^ 00. Then, the first three 
asymptotic moments are given by the following theorem: 

Theorem 1: Assuming that the matrices D and D satisfy the conditions stated in lAll and IA2I then the following 
convergences hold true: 

1) First asymptotic moment [12], [13]: 

and Ef^V- ^0, 



p \P0 / P 

2) Second asymptotic moment [12], [13]: 



VIk = 0{l) and K¥.[^— - n\ 
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3) Third asymptotic moment: 



0. 



VK = 0{l) and 

The two first items of the theorem are proved in [13] (beware that the notations used in this article are the same 
as those in [4] and slightly differ from those used in [13]). Proof of the third item of the theorem is postponed to 
the appendix. 

Remark 1: One can note that the third asymptotic moment is of order 0{K~'^). This is in accordance with the 
asymptotic normahty of the SNR, where the third moment of ^/K{i3k ~ ]E(/3i<-)) will eventually vanish, as this 
quantity becomes closer to a Gaussian random variable. However, its value remains significant for small dimension 
systems. 



V. Simulation results 

In our simulations, we consider a MIMO system in the uplink direction. The base station is equipped with N 
receiving antennas and detects the symbols transmitted by a particular user in the presence of K interfering users. 
We assume that the correlation matrix \£' is given by — •y^^a''"-'' with < a < 1. Recall that P is the 

matrix of the interfering users' powers. We set P (up to a permutation of its diagonal elements) to: 

diag([4P 5P]) if K = 2 
diag([P P 2P 4P]) if 

where P is the power of the user of interest. For K = 2^ with 3 < p < 5, we assume that the powers of 
the interfering sources are arranged into five classes as in Table |V] We investigate the impact of the correlation 

TABLE I 

Power classes and relative frequencies 



Class 


1 


2 


3 


4 


5 


Power 


P 


2P 


4P 


BP 


16P 


Relative frequency 


1/8 


1/4 


1/4 


1/8 


1/4 



coefficient a on the accuracy of the asymptotic moments when the input SNR is set to 15dB for N = K (Fig. [T]) 
and N = 2K (Fig. |2]l. In these figures, the relative error on the estimated first three moments ( and 

fi denote respectively the asymptotic and empirical moment ) is depicted with respect to the correlation coefficient 
a. These simulations show that when the number of antennas is small, the asymptotic approximation of the second 
and third moments degrades for large correlation coefficients (a close to one). Despite these discrepancies for a 
close to 1, simulations show that the BER and the outage probability are well approximated even for small system 
dimensions. Indeed, Figure [3] shows the evolution of the empirical BER and the theoretical BER predicted by ([T]l 
versus the input SNR for different values of a, K and N. In Figure HI the saddle point approximate of the outage 
probability given by (|2|i is compared with the empirical one. In both Figures [3] and |4] 2000 channel realizations 



October 16, 2008 



DRAFT 



8 



have been considered, and in Fig. |4] the input SNR has been set to 15 dB. These figures show that even for small 
system dimensions, the BER is well approximated for a wide range of SNR values. The outage probability is also 
well approximated except for small values of the SNR threshold that are likely to be in the tail of the asymptotic 
distribution. 

Appendix 1 
Proof of Theorem[T] 

In the sequel, we shall heavily rely on the results and techniques developed in [4]. In the sequel, D and D are 
respectively N x N and K x K diagonal matrices which satisfy lAll and IA2i Z is a iV x iiT matrix whose entries 
are i.i.d. standard complex Gaussian, ^ is a N x K matrix defined by: 

X = D^ZD^ . 

We shall often write X = [xi, • • • , x^] where the x^'s are X's columns. We recall hereafter the mathematical tools 
that will be of constant use in the sequel. 

A. Notations 

Define the resolvant matrix H by: 

H = (^-^D^ZDZ*D5 +I^r^ =(J^X.X.*+In 

We introduce the following intermediate quantities: 

(3{t) = — Tr(DH), a{t) ^ — Tr(DEH) and (3 - a . 

K K 

Matrix R(i) = diag (fi, • • • , tk) is a K x K diagonal matrix defined by: 

R(t) = (l + to(i)DA') ^ . 
Let a = -^Tr(DR). Then, matrix R(t) = diag (ri, • • • , rjy) is a x matrix defined by: 

R(t) = {l + ta{t)T>y^ . 




B. Mathematical Tools 

The results below, of constant use in the proof of Theorem [T] can be found in [4]. 
1 ) Differentiation formulas : 



'lit - -^[X-Hl„ff,.^-i:[.-H],i,„, (3, 



dHpq t , t 



= - -^[HX],,^,, = --[Hx,]^i7., (4) 
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2) Integration by parts formula for Gaussian functionals: : Let $ be a complex function polynomially 
bounded together with its derivatives, then: 

'd $(X) 



E[Xy$(X)] 



(5) 



3) Poincare-Nash inequality: Let X and $ be as above, then: 



5$(X) 



9$(X) 



dX,,. 



(6) 



TV K 

Var($(X)) < XI XI ^''^J* 
1=1 j=i 

4) Deterministic approximations and various estimations: 

Proposition 2: Let (A^) and (B^) be two sequences of respectively N x N and K x K diagonal deterministic 
matrices whose spectral norm are uniformly bounded in K, then the following hold true: 



1 



-Tr(AR) 



1 



K ' ' K ^ ' ^ 
Proposition 3: Let {Ak), (Bk) and (Ca-) be three sequences of NxN, KxK and NxN diagonal deterministic 



4Tr(BR) = 4Tr(BT) + 0{K-''). 
K K 



K \ K 



matrices whose spectral norm are uniformly bounded in K. Consider the following functions: 

a>(X) = iTr (aH^^) , .I.(X) 
Then, 

1) the following estimations hold true: 

var$(X), var«'(X), var(/3) and var f -^TrAHCH j are 0{K-^) 



2) the following approximations hold true: 
E [$(X)] 



— Tr f DTb) — Tr (ADT) + 0{K-^), 
K \ J K 



(7) 



E [^'(X)] = , [^^^ (dTb) Tr(AD^T2) - -^Tr (d^T^b) Tr(ADT) j + 0{K-^), (8) 



1 — t'^^^ 

E— TrfAHDHl = \ — — 

K ^ ^ 1 - t^-i^ K 



Tr(ADT2) + 0{K- 



(9) 



Proofs of Propositions |2] and [3] are essentially provided in [4]. In the same vein, the following proposition will be 
needed. 

Proposition 4: Let (A^), (Bj<-) and (C^) be three sequences of NxN, KxK and NxN diagonal deterministic 
matrices whose spectral norm are uniformly bounded in K. Consider the following function: 



^(X) = — Tr 



CHAHAH 



XBX* 

K 



Then var(^(X) = 0{K-^) and var (-^TrAHAHAH) = 0{K-^) . 

Proof of Proposition |4] is essentially the same as the proof of Proposition [3} 1). It is provided for completeness and 
postponed to appendix HIl 
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C. End of proof of Theorem [7] 

We are now in position to complete the proof of Theorem [T] Using the notations of [4], the SNR writes: 

/3k = ^z*D^H(t)D^z, 

where t — j;- Hence, the third moment is given by: 

E^PK-^Pxf = ■^^^E(z*D^HD^z-ETrDH)^ 

- ^^^°^^E (z*D^HD^z - TrDH + TrDH - ETrDH)^ 

E (z*D3HD^z - Ti-Dh)^ + 3E (z*D^HD^z - TiDh)^ (TrDH - ETrDH) 



iC3 

+3E (z*D5HD5z - TrDH^ (TrDH - ETrDH)^ + E (TrDH - ETrDH) 

E (z*D5HD^z - TiDh)^ + 3E (z*D^HD^z - TiDh)^ (TrDH - ETrDH) 

(10) 



f E (TrDH - ETrDH) 



In order to deal with the first term of the right-hand side of ( fTOl i. notice that if M is a deterministic matrix and x 
is a standard Gaussian vector, then: 

E(x*Mx-TrM)^ = Ty{M^)E {\xi\'^ - if 

(such an identity can be easily proved by considering the spectral decomposition of M). Hence, 

E(z*D^HD^z-TrDHy'' = ETr (DH)^ E ([Znp - l)^ 

= 2ETr (DHDHDH) . 

The second term of the right-hand side of ( fTOl ) is uniformly bounded in K. Indeed: 

3E (^z*D^HD^z- Tr(DH))^ = 3E (jZnp - l)^ TrDHDH (TrDH - ETrDH) , 



< 3 Vvar (TrDHDH) ^var (TrDH) 

which is 0(1) according to Proposition [3] It remains to deal with E (TrDH — ETrDH)^, which can be proved to 
be uniformly bounded in K using concentration results for the spectral measure of random matrices [21] (see also 
[15, eq.(86)-(87)], where details are provided). Consequently, we end up with the following approximation: 

K^¥.{(3k - ^Pxf = -^^^EdZiip - l)^ETrDHDHDH + (ivT-i) 

which is deterministic but still depends on the distribution of the entries via the expectation operator E. The rest of 
the proof is devoted to provide a deterministic approximation of ETr (DHDHDH) depending on 7, 7, T and T. 
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Note that H = I - HXX*, thus: 
fHDHDH 



pp 



[HDHD]^p-i 



HDHDH 



XX* 

K 



pp 



[HDHD]pp [HDHDHx,]^ Xp 



(11) 



Let us deal with the second term of (fTTT l. We have: 



1 1 ^ 

E- [HDHDHx,]^ ^^J^T.^ ([HDHDH]^, X^.X, 



k=l 



Using the integration by part formula (|5]l, we get: 

N N 

E [HDHDHxj]p = ^ dkdj5(p - fc)E [HDHDH]pj^ + ^ dudj^ 



k=l 



k=l 



N 



£.m=l 



kj 



N 



dp J,E [HDHDH]pp - - ^ dkdjd^deE X„ [Hx,]^ 



A;,£,m— 1 



^ dkdjdmde^ [XpjHpi [Hxj]^ HkmHmk] 

k,e,m=l 

t ^ 



dpdjE [HDHDH]pp djE [Hxj]^ Xp^-Tr (DHDHDH 



K ^ 



K 

[HDHxj]p^Tr (DHDH) 



- —Le 

K ^ 



[HDHDHxj] XpjTr(DH) 



Substituting in the last term -^TrDH =/? +« where /3= /3 — a, we get: 



E [HDHDHxj]p Xpj = dp J^E [HDHDH]pp - ^dj^ 



[Hxj] XpjTr (DHDHDH) 



-— fl,E 
if ^ 



[HDHxj]p XpjTr (DHDH) 



[HDHDHx,]pXp,/3 



-tdjE 



[HDHDHx,]^Xp, 



Therefore, we have: 



(l + tadj ] E 



[HDHDHx,]pXp, 



iXp, = dpd>[HDHDH]p^--E [Hx,]p Xp, J, Tr [DHDHDH] 



K 



-dM 



[HDHxj]p XpjTi [DHDH] - tdjE 



[HDHDHx 



■ilp^pi P 



October 16, 2008 



12 



Multiplying the right hand and the left hand sides by r 
E[HDHDHxj-]p3;~ = fjdpJjE[HDHDH]pp 



'■> l+tadj 
t ^ 

— r 
K 



, we get: 



[HDHxj] XpjTr[DHDH] 



— tdjfjM 



[HDHDHx,] Xp, /? 



(12) 



Plugging (T% into ( fTTT l. we obtain: 



E [HDHDH]^^ = E [HDHD]^^ - ^ -P,dpd,E [HDHDH]^^ + ]f2 E ^^I^ [Hxj]p Xp.d.Tr [DHDHDH 



J2 djrjE [HDHx,]p Xp,,Tr [DHDH] + -| ^ d.f.E [HDHDHx,]^ Xp,, P, 



E [HDHD]pp - tadpE [HDHDH]pp + -|jETr(DHDHDH) 



+ ^^Tt [DHDH] 



HDHXDRX* 



+ -E/3 



HXRDX* 
HDHDHXDRX* 



pp 



Hence, 

(1 + tadp)E [HDHDH 



- E[HDHD]pp + —ETi [DHDHDH] 



HXRDX* 



pp 



+ — rETr [DHDH] 



HDHXDRX* 



Multiplying the left and right hand sides by r,p = -j^f^^, we get: 

E[HDHDH]pp = TpE [HDHD]pp + -^rpETV [DHDHDH] 



t o 

+ -E/3 

pp I\ 



HXRDX* 



HDHDHXDRX 



pp 



rpETr [DHDH] 



HDHXDRX* 



pp 



K 



rpE/3 



HDHDHXDRX* 



(13) 



pp 



Multiplying by dp, summing over p and dividing by K, we obtain: 
E—Tr [DHDHDH] = E— ^ rfp [HDHDH]^^ , 



= 1 



^ K 2 

— V TpdpE [HDHD]„„ + — E Tr (DHDHDH) Tr f DRHXRDX 



p[ ■ P"P ^ 'pp ' J{3 

m (DHDH) Tr (dRHDHXDRX*) 



p=i 



A 



^E /3Tr( 



X1+X2+X3+ X4, 



DRHDHDHXDRX' 



(14) 
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where: 



Xi 



X2 



X3 



X4 



— ETr (DRHDHD) 

K ^ ' 



ETr (DHDHDH) — Tr DRH 



K 



K 

e_ 

K 

— E (3 Tr DRHDHDH 



ETr (DHDH) — Tr DRHDH 



XDRX* 

K 

XDRX* 



K 



XDRX* 

K 



According to Proposition [3 var-^Tr (^DRHDHDH xdRX* j of order OiR-"^). Similarly, var(/3) = 0{K-'^). 
Hence, using Cauchy-Schwartz inequality, we get the estimation Xi = 0{K^^). It remains to work out the 
expressions involved in xi, X2 and xz by removing the terms with expectation and replacing them with deterministic 
equivalents. 

= 0(i^-2) by Proposition [3] and var(iTrDHDHDH) = 0{K-'^) by 



rXDRX 



Since var^^Tr (^DRH^— ^ 



Proposition m we have: 

X2 = 



11 

—ETr (DHDHDH) E — Tr 

K ^ ^ \ K 



DRH- 



XDRX* 



K 



(a) 



1 



K 



K 



ETr (DHDHDH) — Tr DTDR — Tr (DRDT) + 0{K 



K 



= —ETr (DHDHDH) 77 + 0(i^-2) 
K 



(15) 



where (a) follows from Proposition |3]-2) and (b), from Proposition |2l Similar arguments yield: 

XDRX* 



X3 



i2 I 
—ETr (DHDH) E — Tr 
K ^ ' \ K 



DRHDH- 



K 



0(K- 



(1 - ^2^7)2 

^^7 

(1 - t^iiY 



^Tr (dTDr) -^Tr (DRD^T^) - ^Tr (d^T^Dr) -^Tr(DRDT) 



K 

1 
K 



and 



Tr(D''T 



Xi 



K 



Tr(D'^T'^) 



OiK- 



(16) 



1 1 

1 - i2^7 K 
1 1 
1 - i2^7 K 



Tr(D2RDT^) +e'(X"^) 



Tr(D^T^) +C'(X-2)_ 



(17) 



Plugging ( fT6b . ( fTSl ) and ( fTTI ) into ( fT4b . we obtain: 

— ETr(DHDHDH) = — . n ^, 

K ^ ' i^(l-i2^7)3 /S:(l-t2^7)3 



ETr(DHDHDH) = — ^— — -^TrD^T^ — -^TrT^D^ + OiK-"^). 
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Hence, 



Po Po 



K (p2 _ 



2p3 



TrD^T^ - ^TrD^T^ 



TrD^T^ _ ^TrD^T^ 



Zn\'-lY + 0{^], 



1 



The fact that 



2p' 



Proof of Theorem [T] is completed. 



is of order 0(1) is straightforward and its proof is omitted. 



Appendix II 
Proof of Proposition|4] 

The proof mainly relies on Poincare-Nash inequality. Using the Poincare-Nash inequality, we have: 



N K 

var((^(X)) < ^^rf.Jjl 

i=i j=i 



dip 



dX, 



2 N K 



1=1 j=l 



dip 



dX,, 



We only deal with the first term of the last inequality (the second term can be handled similarly). We have 
p>{%) = Y.p,r.s,t=i Y.u=i CppHprArrHrsAssHstXtuBuuX*^. After Straightforward calculations using the dif- 
ferentiation formula (|3]l, we get that: 

d(p 



dX,,. 



where: 



— - [X*HAHAHXBX*CH1 , 



(2) 



[X*HXBX*CHAHAH],, , 



(4) 



[X*HAHXBX*CHAH1 

-i^[BX*CHAHAH]^,. 



Hence. 







N K 

7 = 1 1 = 1 


dp 
dX^j 



and 



< -^ETr ( DHCXBX* HAH AHXDX* HAH AHXBX* CH 

4 



ETr (^DHAHCXBX*HAHXDX*HAHXBX*CHAH 
ETr (dHAHAHCXBX*HXDX*HXBX*CHAHAH 
ETr ( DHAHAHCXBDBX*CHAHAH 
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We only prove that the first term of the right hand side is of order K ^; the other terms being handled similarly. 
Using Cauchy-Schwartz inequality, we get: 



N K 



J-2J llTJI|2||r'||2 



4^^d,d,E|(/.i/ < ^'"ax^ll ll^-ll ETr((HA)'HXDX*H(AHf (XBX*f J , 

< -r^dra^.^'U.fWCf (ETr (HA)^ HXDX*H (AH)^ (HA)^ HXDX*H (AH) 
X (ETr(XBX*)'')' 



2\ 2 



< |!d_||H||«||C||^||Ar^El(^^°^* 




K I \l K \ K 

where the first inequality follows by using the fact that |TrAB| < ||B||Tr(A), A being hermitian non-negative 



matrix and the second follows by applyig twice Cauchy-Schwartz inequalities: Tr (AB) < \/Tt (AA*)-\/Tr (BB*) 
and EXY < VeX^VeY^. We end up the proof of the first statement by using the fact that -^E [-^Tr (-^XBa-X*)"] 
is uniformly bounded in K whenever B^' is a sequence of diagonal matrices with uniformly bounded spectral norm 
and n is a given integer 

The second statement follows from the resolvent identity: 

—TrAHAHAH = —TrAHAHA — TrAHAHAHXX*. 

K K K 

According to the first part of the proposition, 

var ^-^ TrAHAHAHXX*^ = OiR-"^) . 

Now, TrAHAHA = TrA^HAH and var;iTrA2HAH 0(/\-2) Proposition [3}l). Hence, applying 
inequality var(X + F) < ^&y{X) + var(y) + 2-\/var(X)var(y) yields the desired result. Proof of Proposition |4] 
is completed. 

References 

[1] I. E. Telatar, "Capacity of Multi-antenna Gaussian Channels," European Transactions on Telecommunications and Related Technologies, 

vol. 10, no. 6, pp. 585-596. November 1999. 
[2] B.M. Hochwald, T.L. Marzetta, and V. Tarokh, "Multiple-Antenna Channel Hardening and Its Implications for Rate Feedback and 

Scheduling," IEEE Transactions on Information Theory, vol. 50, no. 9, pp. 1893-1909, September 2004. 
[3] R. Narasimhan, "Transmit Antenna Selection Based on Outage Probability for Correlated MIMO Multiple Access Channels," IEEE 

Transactions on Wireless Communications, vol. 5, no. 10, pp. 2945-2955, October 2006. 
[4] W. Hachem, O. Khorunzhiy. Ph. Loubaton, J. Najim, and L. Pastur, "A New Approach for Capacity Analysis of Large Dimensional 

Multi-Antenna Channels," IEEE Inf. Theory, September 2008. 
[5] G. Taricco, "Asymptotic Mutual Information Statistics of Seperately-Correlated Rician Fading MIMO Channels," IEEE Transactions on 

Information Theory, 2008. 

[6] A.L. Moustakas and S.H. Simon, "On the Outage Capacity of Correlated Multiple-Path MIMO Channels," IEEE Transactions on 

Information Theory, vol. 53, no. 11, November 2007. 
[7] W. Hachem, P. Loubaton, and J. Najim, "A CLT for Information Theoretic Statistics of Gram Random Matrices With a given Variance 

Profile," to be published in Annals of Applied Probability, 2008. 



October 16, 2008 



DRAFT 



16 



[8] Y.C. Ko, M.S. Alouini, and M.K. Simon, "Outage Probability of Diversity Systems over Generalized Fading Channels," IEEE Transactions 

on Communications, vol. 48, no. 11, pp. 1783-1787, November 2000. 
[9] G.L. Stiiber, Principles of Mobile Communications, Norwell, MA:Kluwer, 1996. 
[10] S. Jin, M.R. McKay, X. Gao, and I.B. CoUings, "Asymptotic SER and Outage Probability of MIMO MRC in Correlated Fading," IEEE 

Signal Processing Letters, vol. 14, no. 1, 2007. 
[11] J.W. Stokes and J.A. Ritcey, "A General Method for Evaluating Outage Probabilities Using Pad^ Approximations," Globecom, vol. 3, no. 
8-12, pp. 1485-1490, 1998. 

[12] A. Kammoun, M. Kharouf, W. Hachem, and J. Najim, "A central limit theorem for the snr at the wiener filter output for large dimensional 

signals," IEEE Workshop on Signal Processing Advances in Wireless Communications, 2008. 
[13] A. Kammoun, M. Kharouf, W. Hachem, and J. Najim, "A Central Limit Theorem for the SINR at the LMMSE Estimator Output for 

Large Dimensional Signals," Submitted to IEEE Transactions on Information Theory, 2008. 
[14] D. Guo, S. Verdu, and L.K. Rasmussen, "Asymptotic Normahty of Linear Multiuser Receiver Outputs," IEEE Trans. Info. Theory, vol. 

48, no. 12, pp. 3080-3095, December 2002. 
[15] P. Li, D. Paul, R. Narasimhan, and J. Cioffi, "On the Distribution of SINR for the MMSE MIMO Receiver and Performance Analysis," 

IEEE Inf Theory, vol. 52, no. 1, pp. 271-286, January 2006. 
[16] J. Dumont, W. Hachem, S. Lasaulce, P. Loubaton, and J. Najim, "On the capacity achieving covariance matrix for rician mimo channels: 

an asymptotic approach," submitted to IEEE Inf. Th., 2008. 
[17] D. N. C Tse and O. Zeitouni, "Linear Multiuser Receivers in Random Environments," IEEE Trans. Info. Theory, vol. 46, no. 1, pp. 

171-188, January 2000. 

[18] H. V. Poor and S. Verdu, "ProbabiUty of error in MMSE multiuser detection," IEEE Trans. Info. Theory, vol. 43, no. 3, pp. 858-871, 
May 1997. 

[19] M. K. Simon and M. S. alouini. Digital Communication Over Fading Channels, Wiley Series in Telecommunications and Signal Processing, 

second edition, 2005. 

[20] R. W. Butler and A. T. A. Wood, "Saddlepoint approximation for moment generating functions of truncated random variables," The Annals 

of Statistics, vol. 32, no. 6, pp. 2712-2730, 2004. 
[21] A. GuiUonnet and O. Zeitouni, "Concentration of the Spectral Measitfe for Large Matrices," Electronic Communications in Probability, 

pp. 119-136, 2000. 



October 16, 2008 



DRAFT 



Correlation coefficient 
(a) First moment of the SNR 



25 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Correlation coefficient 



(b) Second moment of the SNR 



60 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Correlation coefficient 



(c) Third moment of the SNR 
Fig. 1. Absolute value of the relative error when N = K 
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Fig. 2. Absolute value of the relative error when = 2K 
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Fig. 4. Outage Probability vs SNR threshold 
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